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Abstract 

The elemental compositions of hot Jupiters are informative relics of planet formation that can help 
us answer long-standing questions regarding the origin and formation of giant planets. Here, I present 
the main conclusions from a comprehensive atmospheric retrieval survey of eight hot Jupiters with 
detectable molecular absorption in their near-infrared transmission spectra. I analyze the eight trans¬ 
mission spectra using the newly-developed, self-consistent atmospheric retrieval framework, SCAR¬ 
LET. Unlike previous methods, SCARLET combines the physical and chemical consistency of complex 
atmospheric models with the statistical treatment of observational uncertainties known from atmo¬ 
spheric retrieval techniques. I find that all eight hot Jupiters consistently require carbon-to-oxygen 
ratios (C/O) below ~ 0.9. The finding of C/O < 0.9 is highly robust for HD 209458b, WASP-12b, 
WASP-19b, HAT-P-lb, and XO-lb. For HD 189733b, WASP-17b, and WASP-43b, I find that the 
published WFCS transmission spectra favor C/O < 0.9 at greater than 95% confidence. I further 
show that the water abundances on all eight hot Jupiters are consistent with solar composition. The 
relatively small depth of the detected water absorption features is due to the presence of clouds, not 
due to a low water abundance as previously suggested for HD209458b. The presence of a thick cloud 
deck is inferred for HD 209458b and WASP-12b. HD 189733b may host a similar cloud deck, rather 
than the previously suggested Rayleigh hazes, if star spots affect the observed spectrum. The ap¬ 
proach taken in SCARLET can be regarded as a new pathway to interpreting spectral observations 
of planetary atmospheres. In this work, including our prior knowledge of H-C-N-0 chemistry enables 
me to constrain the C/O ratio without detecting a single carbon-bearing molecule. 


1. INTRODUCTION 

Spectroscopy observations of hot Jupiters provide an 
unprecedented opportunity to constrain the metallici- 
ties and carbon-to-oxygen ratios (C/O) of giant planets 
which can provide crucial hints on the formation of gi¬ 
ant planets including the gas and ice giants in our solar 
system. If giant planets predominately form in a quick 
one^step process through gravitational instability (Boss 
1997), their atmospheres should have elemen tal ao um 
da nces resem bling the ones of their host star (Helled & 
Schubert 2009). If, on the other hand, gian t planets form 
throu gh multi-step core accretion process (Pollack et al. 
1996), a wide range of elemental compositions with dis¬ 
tinct deviations fr om the host star’s elemental composi- 
tion are expected dOberg et al. 112011 Helling et al. 112014 
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with distance from the star because of the different con ¬ 
d ensa tion temperatures of H 2 O and CO (Oberg et al 
2011). Detailed disk evolution models suggest that the 


Madhusudhan et al.|2014a||Ali-L)ib et al.|2014[ IVlordasini 

et al., 2015). Within the core accretion model, refrac¬ 
tory and ice-forming elements first segregate from the gas 
in the protoplanetary disk through condensation. They 
then evolve separately from the gas for millions of years, 
and finally, some of the gas and solids are recombined in 
the process of planet formation. In this multi-step pro¬ 
cess, it would be a strange coincidence for the final planet 
to form with elemental compositions identical to the one 
of the host star. Instead, we expect the elemental abun¬ 
dances of giant exoplanets to be insightful relics of this 
planet formation process. 

The carbon-to-oxygen ratio (C/O), in particular, has 
been proposed to be a good tracer of planet formation. 
Protoplanetary disk models suggest that the C/O ra¬ 
tio of the gas and solids in protoplanetary disks varies 
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C/O of the gas beyond the water iceline should tran¬ 
sition to C /O -> 1 on time scales of ^ 3Myrs (Helling 
et al.|2 014|). As a result, the gas envelopes of giant plan¬ 
ets” formed 


beyond the iceline may be carbon enriched 
with the exact C to O ratio depending on the time and 
the location of the runaway gas accretion. If, on the 
other hand, gas giant atmospheres are heavily polluted 
by the accretion of oxygen enriched solids, the C/O ratio 
of giant planet envelopes may substantially deviate from 
the C/O ratio in the gas disk. Van Boekel et al. (2015) 
model the formation and evolution of giant exoplanets 
from the formation of planetesimals to today’s evolved gi¬ 
ant planets, keeping track of the amounts of gaseous, icy, 
and rocky material accreted as the planets forms and mi¬ 
grates within the disk. They conclude that giant planet 
formation beyond the water ice line should always lead 
to C/O < 1, while giant planet formation within the ice- 
line can lead to either carbon-rich chemistry (C/O > 1) 
or oxygen-rich chemistry (C/O < 1) depending on the 
carbon abundance of refractory materials. The uncer¬ 
tainties in planet formation and disk evolution models 
are substantial, however. Observational constraints on 
the elemental abundances of giant planets are greatly 
needed. 

The challenge in obtaining reliable observational con¬ 
straints on the elemental abundances lies not only in ob¬ 
taining sufficiently precise data, but also in the difficulty 
associated with inferring elemental abundances represen¬ 
tative of the planet’s bulk gas envelope from spectro¬ 
scopic observations of the planet’s photosphere. Despite 
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exquisite observations, no representative values for the 
C/O ratios of the solar system gas giants, no representa¬ 
tive values for the C/O ratio in the Solar System giants 
have been determined to date, because oxygen is con¬ 
densed out in the form of water in deep layers, largely 
inaccessible to remote sensing. More accurate oxygen 
abundances and C/O ratios for the solar system gas gi¬ 
ants would provide an important new constraint on our 
models of gas giant planet formation. One of the main 
science objectives of the Juno mission — expected to 
arrive in 2016 — is to provide the first representative 
C/O measurement for Jupiter’s gas envelope through mi¬ 
crowave observations from a low Jupiter orbit. 

Observations of hot Jupiters and directly imaged young 
giant planets present a new opportunity to determine the 
elemental abundances for a large sample of giant planets 
using infrared spectroscopy. The temperatures in their 
atmospheres are sufficiently high that all of the dominant 
carbon or oxygen bearing molecular species are accessible 
to IR remote sensing. However, the spectral signatures 
of exoplanets are governed not only by the planet’s ele¬ 
mental composition, but also its radiation environment 
and a wide range of weakly understood chemical and dy¬ 
namical processes in their atmosphere. 

Here, I introduce a new “Self-Consistent Atmospheric 
Retrieval framework for ExoplaneTs” (SCARLET) which 
aims to provide direct constraints on elemental abun¬ 
dances while accounting for the uncertainties resulting 
from our limited understanding of the chemical, dynam¬ 
ical and cloud formation processes in the atmospheres of 
these planets. Unlike previous methods, the new frame¬ 
work is able to determine the full range of self-consistent 
scenarios for a given planet by combining the obser¬ 
vational data and our prior knowledge of atmospheric 
chemistry and physics in a statistically robust Bayesian 
analysis. One key concept is that for a given metallic- 
ity and C/O ratio there is not a single model transmis¬ 
sion spectrum to be compared to the data, but instead 
a range of models corresponding to the uncertainties in¬ 
troduced by our limited understanding of the relevant at¬ 
mospheric processes. SCARLET accounts for this model 
uncertainty by marginalizing over a wide range of “at¬ 
mospheric process parameters” whose prior distributions 
quantify our limited knowledge of vertical mixing, chem¬ 
istry, and cloud formation in the atmospheres of hot ex¬ 
oplanets. 

In this first work, I apply the SCARLET frame¬ 
work to eight published hot Jupiter transmission spec¬ 
tra with detectable near-IR molecular absorption fea¬ 
tures (HD 209458b, WASP-19b, HAT-P-lb, XO-lb, 
HD 189733b, WASP-12b, WASP-17b, and WASP-43b). 
In this study, we focus on HST WFCS observations be¬ 
cause they provide substantial simultaneous spectral cov¬ 
erage (1.1 — 1.7/im), a nd have been shown t o be ex¬ 
tremely repeatable (e.g., Kreidberg et al. 20143), leaving 
little doubt about the reliability of the water detections. 
Where available we include HST STIS and Spitzer ob¬ 
servations in the analysis — note, however, that the C/O 
constraints presented in this work rely only on the HST 
WFCS observations. Our primary conclusions do not 
depend on comparisons of observations taken at differ¬ 
ent times, and the repeatability of WFCS observations 
strongly indicates that the results are not compromised 
by systematic effects that may have affected early C/O 



Fig. 1.— Strength and weakness of self-consistent forward models 
and atmospheric retrieval techniques for the interpretation of ex¬ 
oplanet spectra. Strengths are marked in green. Weakness are 
marked in red. The new “SCARLET” approach combines the 
strengths of self-consistent forward modeling and atmospheric re¬ 
trieval techniques and overcomes their weaknesses. 


measurement s based on broadband Spitzer IRAC data 
points (e.g., iMadhusudh an et ah]|2011 Crossfield et al. 
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In this article, I first motivate the development of the 
new SCARLET framework by reviewing the strengths 
and weaknesses of previous modeling approaches for ana¬ 
lyzing planetary atmospheres (Section 2). Section 3 then 
describes the details of the SCARLET framework. Sec¬ 
tion 4 describes the results from the retrieval study of 
four hot Jupiters. Section 5 presents the conclusions from 
this work and discusses the implications for hot Jupiter 
formation and migration. 


2. BACKGROUND: INTERPRETING PLANETARY 
SPECTRA 

The interpretation of spectral observations of solar 
system atmospheres, exoplanets, and brown dwarfs has 
historically been performed using two distinct mod¬ 
eling strategies: “Atmosphere forward modeling” is a 
theory-driven approach in which the governing equa¬ 
tions of chemistry and physics are solved to derive a 
small number of self-consistent scenarios for the atmo¬ 
sphere. “Atmospheric retrieval”, on the other hand, is 
an observation-driven approach in which statistically ro¬ 
bust constraints on the atmospheric state molecular com¬ 
position and temperature structure are retrieved by re¬ 
peatedly comparing observations to model spectra of pa¬ 
rameterized scenarios without modeling the chemical and 
physical process in the atmosphere. Both approaches 
have strengths and severe limitations (Figure IT]). The 
goal of this section is to provide a concise overview of “for¬ 
ward models” and “atmospheric retrieval techniques” to 
motivate the development of the new SCARLET frame¬ 
work presented in this work. 


2.1. Complex Atmospheric “Forward” Models 

The basic idea of atmosphere forward models is to it¬ 
erate the relevant atmospheric processes including atmo¬ 
spheric chemistry, radiative-convective heat transport, 
and/or cloud format ion until a converged self-consistent 
solution is achieved ( Burrows et al.|[l997 Seager & Sas- 
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(e.g., 

Zahnle et al. 2009 Moses et al. 201 la| Hu et al. 

2012 

2009 

, the atmospheric dynamics (e.g., Showman et al. 

, the thermal structure (e.g., McKay et al. 1989; 

Burrows et al. 1997 Fortney et al. 2011} and/or the for- 

mation of clouds (e.g., IVfarley et al. 2002 Lavvas et al. 

2010 

Morley et al. 2013). The main strength of these at- 


insights into the atmospheric processes at play. Match¬ 
ing the observations generally requires modeling all rele¬ 
vant physical processes. If solved iteratively, atmospheric 
scenarios from atmosphere forward models are generally 
physically self-consistent within the modeled physics. 

Severe limitations of atmosphere “forward” modeling 
are, however, that they inherently assume that all phys¬ 
ical and chemical processes are known without uncer¬ 
tainty (Figure]!]). Mismatch between data and models is 
often reduced by manually tweaking individual model pa¬ 
rameters until a sufficient match to the data is achieved. 
Manual tweaking, however, prevents a statistically robust 
treatment of the observational uncertainty. Individual 
scenarios matching the data are identified without un¬ 
derstanding what the degeneracies are and whether the 
found scenario is the only plausible scenario. 



Posterior constraints 

"Formation" Parameters: 

• Metallicity 

• C/O 

• Residual heat, T int 

"Atmospheric Process" 

Parameters 

• Eddy diffusion, K zz 

• Bond albedo, A B 

• Day-night heat transport 

• Cloud properties 


Fig. 2. — Schematic view of the SCARLET framework. 


et al. 2013a)). Equally important, atmospheric retrieval 
techniques for exoplanets require an oversimplified de¬ 
scription of atmosphere in order to reduce the number of 
free parameters. Well-mixed atmospheres composed of 
only a small number of pre-selected molecules are con¬ 
sidered, even though the mixing ratios real istic atmo 


spheres can be strongly altitude dependent (Moses et al. 
2013b a). The presented shortcoming of atmospheric re- 
trieval techniques suggest the potential for an improved 
way to analyze spectroscopic observations of exoplanets 
— one that combines the statistical robustness of at¬ 
mospheric retrieval techniques with the physical insights 
provided by atmosphere forward models. 


2.2. Atmospheric Retrieval 

Atmospheric retrieval methods are different from at¬ 
mosphere “forward” models in that they generally do not 
model any of the chemical and dynamical processes in the 
atmosphere. Instead, they describe the molecular com¬ 
position and temperature profile by free parameters and 
model the planet’s spectrum for comparison to the data. 
Atmospheric techniques are powerful when good data is 
available and have a long history of app lications in mete¬ 
orology (e.g. Smith||1970 Chahine||1974 ) and the explo - 
ration of the solar system planets (e.g., Hanel et al. 2003). 
An excellent overview of the theory of atmospheric re- 
trieval within the ass um ptions of Gaussian uncertainties 
is provided in Rodgers (2000). Techniques that more 
robustly explore the full para meter space of p lan etary 
atmos p heres are described in Madhusud han fc Se ager 
( 2009]) ; Benneke & Seager (2012|); Line et al. (2012). 

A key disadvantages of atmospheric retrieval tech¬ 
niques is that they provide no direct insights into the 
physical and chemical processes at play. They also make 
no use of our prior substantial knowledge of atmospheric 
chemistry and heat transport to better constrain the 
properties of the planet studied. This is particularly criti¬ 
cal for exoplanets because the available data is sparse and 
including our prior knowledge of chemistry or heat trans¬ 
port can substantially reduce the uncertainties in the in¬ 
ferred elemental abundances and internal heat flux. In¬ 
stead, atmospheric retrieval techniques make no attempt 
to rule out unphysical scenarios that grossly contradict 
our basic understanding of atmospheric chemistry and 
radiative-convective heat transport. Combinations of 
molecular gases that are highly reactive with each other 
at given cond itions are not ruled out in regular retrieval 
methods (see[Hu & Seager 20l4|for a discussion). Sim¬ 
ilarly, no consistency check between the parameterized 
temperature-pressure profile and the m olecu lar compo¬ 
sition is performed (Madhusudhan & Seager 2009 Line 


3. METHODOLOGY: THE SCARLET FRAMEWORK 

The main objective of the SCARLET framework is 
to determine the full range of self-consistent scenarios 
for an observed planet by combining the observational 
data and our prior knowledge of atmospheric chemistry 
and physics in a statistically robust Bayesian analysis. 
SCARLET consists of a state-of-the-art self-consistent 
photochemistry model surrounded by the nested sam¬ 
pling algorithm for the retrieval. This architecture en¬ 
ables SCARLET to provide direct insights into the ele¬ 
mental composition of the deep atmosphere and plan¬ 
etary bulk based on spectral observations of the up¬ 
per atmosphere. Compared to previous modeling ap¬ 
proaches, SCARLET combines the statistical robustness 
and exploratory nature of traditional atmospheric re¬ 
trieval methods (Smith et al., 1970; Rodgers et al. 2000; 
Madhusudhan & Seager 2009; Benneke & Seager 2012) 
with the self-consistency and ability to provide physical 
insights of so-called self-consistent atmospheric forward 
models (e.g., Burrows et al. 1997, 2001; Marley et al. 
2002; Fortney et al. 2008, 2011; Moses et al 2013). 

SCARLET is distinct from any previous atmospheric 
retrieval methods that describe the current state of the 
visible atmosphere (molecular abundance and tempera¬ 
ture profiles) by free parameters. The advantage of the 
SCARLET approach are as follows: (1) SCARLET pro¬ 
vides direct insight into the elemental composition of 
the planet’s deep atmosphere, (2) all atmospheric sce¬ 
narios found by SCARLET are inherently self-consistent 
(within the limits of our prior knowledge of atmospheric 
chemistry and physics), (3) SCARLET can provide direct 
constraints on physical processes such as the strength of 
vertical mixing in the studied atmospheres. 

3.1. SCARLET Overview 

Figure [2] illustrates the basic architecture of the SCAR¬ 
LET framework. The input to the framework are the 
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TABLE 1 

“Formation” parameters (top) and “atmospheric process” 

PARAMETERS (BOTTOM) AND THEIR PRIOR RANGES IN THE SCARLET 
RETRIEVAL ANALYSIS 


SCARLET 

Prior range 

Metallicity at 1000 bar 

0.1... 100 

C/O at 1000 bar 

0.0001...10000 

Internal heat 

0...200K 

Eddy diffusion coefficient, K zz 

10 7 ... 10 11 cm 2 s —1 

Bond albedo 

o 

o 

lo 

Dayside-nightside heat redistribution factor 

0.25... 0.5 

Cloud properties 

see Section 3.3.3 


observed planetary spectrum and an explicit, quantita¬ 
tive description of our prior knowledge of atmospheric 
chemistry and physics. The framework then repeat¬ 
edly runs a self-consistent atmospheric chemistry and 
radiative-convective model to systemically explore atmo¬ 
spheric scenarios to derive posterior constraints on the 
underlying elemental composition of the planet’s gas en¬ 
velope and the strength and efficiencies of process at play 
in the atmosphere. 


3.2. Planet Parameterization 

The fundamental idea of SCARLET is to divide all 
unknown properties of the planet into “formation his¬ 
tory” parameters and “atmospheric process” parameters. 
Formation history parameters describe properties of the 
planet that are set by the planet’s formation and evolu¬ 
tion history, such as its elemental composition and resid¬ 
ual accretion heat. Atmospheric process parameters de¬ 
scribe the strengths or efficiencies of ongoing chemical 
and physical processes in the atmosphere. 


3.2.1. “Formation History” Parameters 

One of the main goals of studying exoplanets is under¬ 
standing the formation and evolution history of planets. 
Planet properties that provide hints on the planet’s for¬ 
mation histories are the bulk elemental composition of 
the planet’s gas envelope and the residual interior accre¬ 
tion heat. While spectroscopic observations of exoplan¬ 
ets cannot probe these properties directly, SCARLET 
models the chemistry, radiative transfer, and mass and 
heat exchange between the observable optically-thin at¬ 
mosphere and the planet deep atmosphere in order to 
derive conclusion on the formation history parameters 
based on spectroscopic observations of the planet. 

Elemental Composition — In this work, I parameterize the 
elemental composition of the deep atmosphere by the 
metallicity (M) and carbon-to-oxygen ratio (C/O) at the 
1000 bar level. The metallicity is defined as 


M = 


[X/H] 

[X/H] 

Solar 


(i) 


where H is the abundance of hydrogen by number and 
X is abundance of any element heavier than helium. A 
metallicity of 10 times solar (M = 10) and the fiducial 
C/O = (C/0) Solar = 0.54 means that all elements heav¬ 
ier than Helium are ten times more abundant than they 
are in the Sun. A log-uniform prior is assigned for metal- 
licities between 0.1 and 100 (Table [TJ. The helium to 
hydrogen ratio is assumed to be equalto the solar value. 


The C/O parameter additionally allows for variations 
in relative abundances of carbon and oxygen, while keep¬ 
ing the sum of the mixing ratios of carbon and oxygen 
constant. C/O ratios between 10 -5 (100,000 times more 
C) and 10 5 (100,000 times more O) are explored. A 
custom-made stretching is applied to ensure an effective 
exploration of the parameter space between C/O = 0.1 
and C/O = 10, where most of the changes in chemistry 
take place (see y-axis in Figure [8]). A uniform prior prob¬ 
ability is applied in the stretched parameter space, cor¬ 
rectly reflecting our prior understanding that scenarios 
with extreme C/O ratios are less likely than C/O ratios 
between 0.1 and 10. 

Interior heat — The interior heat is described by the ef¬ 
fective internal temperature, 

T «=(Si; f' < 2) 

where Frotai is the total cooling power of the planet 
interior, 47rRp is the planet’s surface area at the top of 
the atmosphere, and a is the Stefan-Boltzmann constant. 
A uniform prior on the log scale is between 0 K and 
200 K. 


3.2.2. u Atmospheric Process” Parameters 

The atmospheric process parameters describe chemical 
and physical processes in the atmosphere that are weakly 
understood a priori, but highly relevant for the interpre¬ 
tation of the observed spectrum. In this work, the at¬ 
mospheric process parameters include the eddy diffusion 
coefficient (K zz ), the properties of clouds, the dayside- 
nightside heat transport. Atmospheric process param¬ 
eters are essential for constraining the desired elemen¬ 
tal composition because they provide the connection be¬ 
tween the desired elemental compositions and the observ¬ 
able near-infrared spectrum of the studied hot Jupiters 
and enable me to account for our uncertain prior knowl¬ 
edge of atmospheric processes. 

Eddy diffusion — The eddy diffusion coefficient (K zz ) is 
an essential free parameter because the vertical transport 
of atmospheric constituents in hot Jupiter atmosphere is 
weakly understood and it sensitively affects how closely 
the observable upper atmosphere of the planet resembles 
the composition in the deep atmosphere. UV photons 
from the star break up molecular bonds in the upper at¬ 
mosphere, thereby significantly reducing the molecular 
abundances of species such as methane and water va¬ 
por in the upper atmosphere. The level to which their 
abundances are reduced depends sensitively on the ver¬ 
tical transport of atmospheric constituents that contin¬ 
uously replenish these abundances. Deriving constraints 
on the bulk elemental composition of the gas envelope 
from molecular absorption formed in the upper atmo¬ 
sphere, therefore, requires me to account for the uncer¬ 
tainty introduced by the unknown eddy diffusion coeffi¬ 
cient. For hot Jupiters, eddy diffusion coefficients are un¬ 
certain over several orders of magnitude for hot Jupiters 
(K zz = 10 7 — 10 10 cm 2 s -1 ). 

Heat redistribution — The dayside-nightside heat redistri¬ 
bution is parameterized through a factor, /, that relates 
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Fast track 




Fig. 3. — Flow chart of the self-consistent photochemistry- 
thermochemistry model applied in SCARLET. 

the irradiation of the ID model atmosphere to the stel¬ 
lar flux at the planet’s orbital distance from its host star. 
The extreme limits are full heat redistribution, i.e. the 
incoming radiation is uniformly distributed across the 
planet, and no dayside-night side heat transport, i.e. all 
the radiation is deposited on the planet’s dayside. The 
heat redistribution factor ranges from / = 0.25 for ef¬ 
ficient heat redistribution to / = 0.5 for no dayside- 
nightside heat transport. 

Clouds — Two different cloud modeling approaches are 
considered in this work. The “parameterized particle size 
and cloud profile” model captures the wide range of cloud 
effects in exoplanet transmission spectra through four 
parameters describing the effective grain size (r e ff), the 
cloud top pressure (ptop)> the condensate mole fraction 
(g*), and a cloud profile shape factor (H c ). Alternatively, 
a simplified, two-parameter “Gray Clouds + Rayleigh 
hazes” model describes the clouds in the atmospheres by 
specifying the cloud top pressure of a gray cloud deck and 
the additional Rayleigh-like opacity due to small parti¬ 
cles. The details of the both c loud modeling approaches 
are discussed in Section 13.3.31 


is a chemical equilibrium model minimizing Gibb’s free 
energy in each layer of the atmosphere; the high com¬ 
plexity model is a full thermo- and photochemical kinet¬ 
ics and transport model. Both models were specifically 
designed to maximize computational efficiency to enable 
thousand of evaluations necessary for the parameter ex¬ 
ploration in SCARLET. The advantage of the chemical 
equilibrium model is its computational efficiency. Chem¬ 
ical equilibrium abundances can be computed > 10 3 
times faster than performing a full evaluation of the 
the thermo- and photochemical kinetics and transport 
model. Inferring the elemental abundance of the deep 
atmosphere based on the observable upper atmosphere, 
however, requires capturing transport-induced quenching 
in the mid-atmosphere, and photochemistry in the upper 
atmosphere, which can only be captured by a full thermo- 
and photochemical kinetics and transport model. 


Equilibrium chemistry model — The low-complexity equi¬ 
librium chemistry model computes the molecular com¬ 
positions from the elemental abundances by minimiz¬ 
ing Gibb’s free energ y in each layer of the atmosphere 
(White et ah , 1958; Miller-Ricci et al]|2009 Benneke & 


Seager 2013). I consider approximately 200 molecular 


species from the elements H, He, C, N, O, Ne, Na, Mg, 
Al, Si, P, S, Cl, Ar, K, Ca, Ti, Cr, Mn, Fe, and Ni. 
Rapid computation of ^ 10 5 molecular compositions is 
accomplished by interpolating the molecular abundance 
for each species from a precomputed 4-dimensional equi¬ 
librium grid (30x30x60x30) spanning the dimensions at¬ 
mospheric metallicity, C/O ratio, pressure, and temper¬ 
ature. 


Kinetics and transport model — The thermo- and pho¬ 
tochemical kinetics and transport model computes the 
molecular composition of the atmosphere by solving the 
coupled mass-continuity equations as a function of pres¬ 
sure for each molecular species: 


3.3. Self-Consistent Photochemistry-Thermochemistry 

Model 


dnj d$j 
dt dz 


— Pi Li 


( 3 ) 


The objective of the photochemistry-thermochemistry 
model is to compute self-consistent model atmospheres 
and synthetic planetary spectra for any plausible elemen¬ 
tal composition, internal heat flux, stellar irradiation, 
and atmospheric process parameters. Self-consistent 
molecular compositions and temperature structures are 
obtained by iteratively solving atmospheric chemistry, 
radiative-convective heat transport, and hydrostatic 
equilibrium. A unique feature of the model is that, in 
its most sophisticated form, it self-consistently solves 
the atmospheric chemistry via a full thermo- and photo¬ 
chemical kinetics and transport model and the radiative- 
convection heat transport via a high-resolution, line-by- 
line radiative transfer scheme (Figure [3]). This approach 
enables the derivation of self-consistent molecular com¬ 
positions and temperature structures for any plausible 
elemental composition. The model complexity can be 
reduced by approximating the composition through com¬ 
putationally efficient thermal equilibrium calculations. 

3.3.1. Chemical Models 

Two chemical models with different levels of sophisti¬ 
cation are used in this work. The low-complexity model 


where ni is the number density of species i (m -3 ), is 
the vertical transport flux of species species i (m -2 s -1 ), 
Pi is the chemical production rate (m -3 s -1 ), and Li is 
the chemical production rate (m -3 s -1 ), and z is the ver¬ 
tical coordinate in the atmosphere (m). The chemical 
production and loss rates are computed by summing the 
contributions of 1760 chemical reactions given the molec¬ 
ular abundance, pressure, temperature in each layer. The 
vertical transport flux includes eddy diffusion and molec¬ 
ular diffusion 

<!>,: = -K zz N^-DN^f+Dn ( 

oz oz \ 

( 4 ) 

where K zz is the parameterized eddy diffusion coeffi¬ 
cient, N is the total number density, fi = ni/N is the 
mixing ratio, D is the molecular diffusion coefficient, Hq 
is the mean scale height, Hi is the scale height of species 
z, and olt is the thermal diffusion factor. 

The kinetics and transport model, in principle, fully 
captures the three main chemical processes in plane¬ 
tary atmospheres—thermochemical equilibrium in the 
deep atmosphere, transport-induced quenching in the 


1 1 a T dT\ 

lfo + W i + ~Y~dz ) ’ 
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mid-atmosphere, and photochemistry in the upper atmo¬ 
sphere. A state-of-the-art reaction list of 1760 chemical 
reactions for 92 molecular species formed by the elements 


H, C, O, and N is adopted from Moses et al. (2013a) 
Thermochemical equilibrium is recovered at high pres^ 
sure and temperature because the chemical reaction list 
includes a matching reverse reac tion for each chemical 
reaction ( Visscher fc Moses||2011 ). 


For the boundary conditions, 1 assume a zero-flux con¬ 
dition at the top of the atmosphere boundary and im¬ 
pose the chemical equilibrium abundances for the speci¬ 
fied elemental abundance of the planet’s gas envelope at 
the lower boundary. The use of a zero-flux condition at 
the upper boundary assumes that the atmospheric loss is 
negligible; Moses et al. (2011) discuss this assumption in 
detail and show that it has no effect on the atmospheric 
composition in the photosphere. The lower boundary is 
set at 1000 bar where temperatures are above 2500 K 
such that chemical equilibrium mixing ratios will pre¬ 
vail. The purpose of SCARLET framework is to con¬ 
strain the elemental abundances at the lower boundary, 
i.e. the boundary conditions to the kinetics and transport 
model, based on observations of the upper atmosphere. 
At the top of the atmosphere, the computational domain 
reaches to 10“ 10 bar; this is sufficiently high to capture 
the UV dissociation of all major molecular species. The 
effect of the zero-flux condition at the upper bo undary 
neglec ts atmo spheric loss; throughly discussed in Moses] 


et al.| ( 2011 a), showing that neglecting atmospheric loss 
nas no effect on the stratospheric or lower-thermospheric 
results. 

Equations [3] and [4] form a stiff system of 9200 coupled 
differential equations for a 100 layer atmosphere, pre¬ 
senting a substantial computational challenge. The stiff¬ 
ness of the system results from the vastly different reac¬ 
tion times scales of the modeled reactions. To efficiently 
solve the stiff system, I use a quasi-const ant step-size im¬ 
plementation of the Numerical Differentiation Equation 
(NDFs). NDFs are found to be more efficient for chem¬ 
ical kinetics problems as compared with inverse Euler 


chemistry codes 

(Allen et al. 1981 

Hu et al. 

2012 y 

loses et al. 2013a|. 


ing the J acobian matrix only when it has substantially 


changed (Brown et al. 1989). Despite these improve¬ 
ments, the use of the full the full kinetics and transport 
model remains the primary bottleneck for SCARLET re¬ 
trieval runs. Evaluation of the kinetics and transport 
model are initialized with a converged solution from iter¬ 
ating the chemical equilibrium and radiative convective 
model (Figure [3]). 

3.3.2. Line-by-Line Radiative-Convective Model 

The purpose of the radiative-convective model is 
to compute temperature-pressure profiles self-consistent 
with the stellar insolatio n and molecular composition of 
the atmosphere (Section |3. 3.1 1). To accomplish this, the 
radiative-convective model iterates the T-p profile until 
the radiative downward flux F ^ due to direct star light 
and infrared reemission is matched by the upward emis¬ 
sion flux in the atmosphere . In the process, atmo¬ 
spheric layers with temperature gradients exceeding the 
adiabatic lapse rate — ^ > T = -M- are declared convec¬ 


tive and adjusted to the adiabatic lapse rate. 

Converging to radiative-convective equilibrium can re¬ 
quire hundreds of radiative transfer calculations of the 
entire atmospheric column. In this work, I present a 
novel numerical technique to speed up repeated line-by- 
line radiative transfer calculations by orders of magni¬ 
tudes (Appendix A). The new numerical technique en¬ 
ables the efficient computation of self-consistent T-p pro¬ 
files for any molecular composition based on line-by- 
line equivalent radiative transfer. Line-by-line radiative 
transfer has the advantage that it accurately describes 
the radiative transport for any atmospheric composition, 
while analytical solution for (semi-)gray opacities can 


d eviate significantly fro m the exact solution (e.g., Par- 


mentier & Guillot 2014) and the assumption of corre¬ 


lated opacities across the entire atmospheric column can 
also become problematic, in particular for atmospheres 
with temperature-pressure profiles near the CH 4 /CO or 
N 2 /NH 3 transition for which the mixing ratios of domi¬ 
nant absorbers can be strongly altitude-dependent. 


3.3.3. Cloud Models 

Clouds play a dominant role in shaping the observable 
transmission spectra of exoplanets, yet the formation and 
chemistry of clouds on exoplanets are poorly understood. 
As a result, clouds present a substantial uncertainty in 
retrieving the desired atmospheric composition of exo¬ 
planets. Given the lack of a reliable predictive model for 
clouds on exoplanets, the approach taken in SCARLET 
is to describe the clouds properties using free parame¬ 
ters with the goal of capturing the full range of plausible 
cloud properties. Marginalizing over the cloud properties 
allows us to determine robust constraints on the planet’s 
gas envelope composition while accounting for our lim¬ 
ited prior understanding of cloud formation. Two cloud 
parameterization are applied in this work: 


“Parameterized Cloud Profile and Particle Size” Model — 
The “parameterized particle size and cloud profile” model 
captures a broad range of cloud effects in exoplanet trans¬ 
mission spectra by exploring different cloud compositions 
and describing the particle size and vertical cloud density 
profile via free parameters in the fit. Capturing the up¬ 
permost structure of the clouds is particularly important 
because transmission spectroscopy probes relatively low 
pressures as starlight travel s on a slant path through the 
planet’s upper atmosphere (|Fortney & Marley 2005) 

The approach taken here is to describe the vertical 
cloud density profile of the uppermost clouds or hazes 
by 


Qc 0 ) = Q* (log P - logptop )^ 0 forptop <P< Phase (5) 

where q c = nc 0 nd/^H 2 is the condensate mole fraction at 
pressure p, He is the cloud profile shape factor, p top is 
the cloud top pressure, Pbase is the pressure at the cloud 
base, and g* is the condensate mole fraction one scale 
height below the cloud top (Figure [4]). The functional 
form of Equation [5] is chosen to resemble cloud profiles 
observed and modeled for the solar system and brown 
dwarfs. The particle size distribution is described by a 
log-normal distribution 


n ( r ) 


\ _ (log a-log r eff ) 2 

- e 2(log crp 

\[FkV log G 


(6) 





































7 



Ptop 
6 ■ Ptop 


-2 


Fig. 4. — Example cloud profiles for the “parameterized cloud 
profile and particle size” model. The parameterization encom¬ 
passes a wide range of cloud profiles including cloud profiles re¬ 
sembling th e ones of typical co ndensation clouds (red, blue curves 
— compare |Ackerman & Marley| {2001 )), cloud profiles that as¬ 
sume a uniform condensate mole fraction (green, dashed line — 
compa re lEtangs et al.l J2 0 08l)h and sharp cloud decks (yellow 
dashed[Benneke <V beager| (|2(J12|). The motivation for parame¬ 
terizing tJie clouds through tree parameters, rather than modeling 
them self-consistently, is that our prior knowledge of the nature 
and formation of clouds on hot Jupiters is generally insufficient to 
reliably capture even the basic trends in the cloud formation. 


where r e ff is the parameterized effective particle size 
and the effective variance is set to cr = 2 . Altogether, pa¬ 
rameterized particle size and cloud profile describes the 
uppermost cloud deck by the four free parameters (q c , 
p c , i7 c , and p top ). Transmission spectra generally con¬ 
vey little information about the cloud base, thus Pbase is 
assumed to be high enough that it has a minimal effect 
on the transmitted starlight. Marginalization over the 
wide range of cloud profiles considered in the “parame¬ 
terized grain size and cloud profile” model ensures that 
retrieved constraints on the atmospheric composition are 
largely independent of assumptions on the unknown ver¬ 
tical transport of cloud particles. 

Particles composed of MgSi0 3 , MgFeSi 04 , and SiC 
are considered for the hot Jupiters studied in this work 
(T eq = 1410 ... 2010). The radiative pr operties are c om¬ 


puted using Mie scattering theory (Hansen & Travis 


1974D and the compl ex refractive indices taken from 


Dorschner et al. (1995). 


Gray Clouds + Rayleigh Hazes — Most published obser¬ 
vations of exoplanets provide little information to obser- 
vationally constrain the detailed cloud profiles or grain 
size. Given sparse data, it can be instructive to investi¬ 
gate the constraints on cloud top pressure and haze opac¬ 
ity from a low-complexity “gray cloud deck + Rayleigh 
hazes” model. The two-parameter “gray cloud deck + 
Rayleigh hazes ”cloud model simultaneously allows for 
the presence of a gray cloud deck as well as “Rayleigh” 
hazes composed of small particles (r p <C A). 

Gray clouds are modeled as a sharp cutoff to the 
planet-grazing starlight below a parameterized pressure 
level in the atmosphere. A sharp cutoff approximates the 
appearance of an upper cloud deck composed on large 
particle or a sudden increase cloud density. Rayleigh 
hazes are assumed to be composed of small particles 
(r p <C A). The opacity of the particles follows k oc 

&o .4 (0-4 /im/A) -4 , where &o .4 is the extinction coefficient 
(cm 2 /g) at 0.4 pm. 


3.3.4. Opacities 

The radiative transfer calculations in this work include 
opacities due to molecular absorption, collision-ind uced 
broad ening from H 2 /H 2 and H 2 /He collisions (Borysow 
2002 h , Rayleigh scattering (see|Benneke & Seager 2012 p, 
and Mie scattering of cloud and haze particles. Molecu- 
lar absorption cross sections are determined directly from 
the molecular line lists provide d in t he hig h-te mperature 
ExoMol database phnnyson & Yurchenko 2012 ]) for CH 4 , 
NH 3 , and TiO, HIlEMF database ([Rothman et al.|2010) 
for H 2 Q, CO, and CO 2 , and the H 1 TRA 1 N database 


(Rothman et al. 2009) for O 2 , O3, OH, C 2 H 2 , C 2 H 4 , 
<J 2 H 6 , H 2 O 2 , ana HO 2 . Absorption by the alkali met¬ 
als (Li, Na, K, Rb, and Cs) is modeled based on the 


line st rengths provided in the VALD database (Piskunov 
et al. |1995) and using the ^-broade ning prescription 


(2003) 


provided in Burrows & Volobuyev 

To speed up the evaluation of a large number of atmo¬ 
spheric models, I precompute the wavelength-dependent 
molecular cross sections for each of the considered molec¬ 
ular species on a temperature and log-pressure grid and 
then interpolate the cross section for the required con¬ 
ditions. Similarly, the scattering properties of Mie scat¬ 
tering of cloud and haze particles are precomputed as a 
function of particle size and wavelength. 

In the upper atmosphere, molecular absorption lines 
become increasingly narrow, requiring a very high spec¬ 
tral resolution to exactly capt ure the shapes o f the thin 
Doppler-broadened lines (Goody & Yung 1995). Instead 


of ensuring that each line shape at low pressure is repre¬ 
sented exactly, I choose an appropriate spectral resolu¬ 
tion for the line-by-line simulation by ensuring that the 
simulated observations are not altered by more than 1 % 
of the observational error-bar when the spectral resolu¬ 
tion is doubled or quadrupled. 


3.3.5. Transmission Spectrum Model 

Once the chemistry-radiative convective model has 
converged to a steady-state solution, I compute model 
transmission spectra and synthetic instrument outputs 
to evaluate the fit between the model and the observa¬ 
tions. The model computes the absorption and scattering 
of stellar light by the planetary atmosphere as the rays 
traverse the day-night terminator region. Extinction due 
to molecular absorption, Rayleigh scattering, and cloud 
particles absorption and scattering is accounted for. In 
the parameterized cloud profile, the optical properties of 
finite-size cloud particles are computed using Mie scat¬ 
tering theory. In the “gray cloud deck + Rayleigh haze” 
model, the gray cloud deck is modeled as a sharp cut¬ 
off of transmission below the parameterized cloud deck 
pressure level. Rayleigh hazes are included as an ad-hoc 
opacity source with a oc A -4 . Finally, the high-resolution 
transmission spectrum is integrated over the instrument 
response function of the individual instrument channels 
for comparison to the astronomical observations. 


3.4. Nested Sampling for Atmospheric Retrieval 

SCARLET employs the multimodal nested sampling 
algorithm, MultiNest, to efficiently explore the multidi¬ 
mensional parameter space and compute the posterior 
distribution of the formation and atmospheric process 
parameters. The mathematical details of the nested sam- 
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Fig. 5.— Joint constraints on the metallicity and C/O ratio in the deep atmospheres of HD 209458b, WASP 19b, WASP-12b, and 
HAT-P-lb. The coloring indicates the normalized posterior probability marginalized over all remaining atmosphere and cloud parameters. 
Black contours mark the 68% (lcr), 95% (2a), and 99.7% (3a) Bayesian credible regions. The HST WFC3 transmission spectra of all 
four planets reveal robustly oxygen-dominated (C/O < 0.9) atmospheres. Constraints on the metallicity, i.e. the overall fraction of species 
heavier than He relative to solar, are weak. Results are sh own for a cloud model in which the particle size and vertical cloud density profile 
are described by free parameters in the fit (Section |3. 3. 3|). The posterior distribution is marginalized over the cloud parameters, the eddy 
diffusion coefficient, the planetary Bond albedo, the dayside-night side heat redistribution. 


pling algorithm are described i n Skilling 2004 Feroz et al 


2009; Benneke & Seager 2013|and its first application to 


atmos pheric retrieval is discussed in |Benneke fe Seager 


(2013). Here, I provide overview of the technique and its 
unique capabilities specifically for atmospheric retrieval. 

Similar to the widely used Markov Chain Monte Carlo 
(MCMC), nested sampling is a Monte Carlo approach to 
efficiently compute the multidimensional posterior dis¬ 
tribution of model parameters in parameter estimation 
problems. The main advantage of nested sampling for 
atmospheric retrieval is, however, that nested sampling 
reliably captures highly no n-Gaussian and multimodal 


posterior distrib utions (e.g. Feroz et al. 2009 Benneke 


& Seager 2013). This is extremely important for at¬ 


mospheric retrieval because strong correlations and de¬ 
generacies between atmospheric parameters are common 
due to the generally sparse data available for exoplanets 
and the complex way in which vital information about 
planet is encoded in the observable planetary spectrum. 
Importantly, nested sampling also provides an excellent 
overview of the goodness-of-fit all across the entire pa¬ 
rameter space, even in regions far away from the maxi¬ 
mum probability (e.g., Figure [8]). This provides an im¬ 
proved ability to assess which regions of the parameters 
space can be excluded at levels of significance higher than 
the commonly stated 95% or 99.7% probability. 

Nested sampling begins by randomly taking a user- 


specified number of “active samples” [N fs 100 ... 10000) 
from the entire prior parameter space. The active sam¬ 
ples are then migrated towards regions of high likelihood 
by repeatedly replacing the lowest-likelihood sample by 
a new sample with a likelihood higher than the sample 
to be rep. In this “outside-in” approach, the complete 
prior parameter space in traversed providing high confi¬ 
dence that the global minimum is identified and that ex¬ 
tremely correlated or banana-shaped posteriors are cap¬ 
tured correctly. Replaced samples are stored such that a 
relatively good description of the probability across the 
entire parameters space can be inferred. In this work, I 
employ simultaneous ellipsoidal nested sampling to effi¬ 
ciently find sample with likelihoods hig her t han the pre¬ 
viously lowest-likelihood sample ( Feroz et al.|[20 09). 


4. RESULTS 

In this section, I present constraints on the atmo¬ 
spheric composition and cloud properties in eight hot 
Jupiters (HD 209458b, WASP-19b, HAT-P-lb, and XO- 
lb, HD 189733b, WASP-12b, WASP-17b, and WASP- 
43b) based on their measured transmission spectra. To 
capture the uncertainty in the atmosphere composition 
introduced by the broad variety of plausible cloud prop¬ 
erties, I present constraints on the composition based on 
the “parameterized particle size and cloud profile” model. 
Cloud-introduced uncertainties on the composition are 
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Fig. 6.— Same as Figure [5] but for the hot Jupiters XO-lb, HD 189733b, WASP-17b, and WASP-43b. The HST WFCS transmission 
spectra of all four planets favor C/O < 0.9 atmosphere at greater 95% confidence. A low-probability tail towards high C/O remains for 
WASP-17b, and WASP-43b containing a few percent of the posterior probability. To illu strate the tail, the probability is plotted over the 
full range from C/O = 10 -5 to C/O = 10 +5 . ( Madhnsndhan et al.|2011 Stevenson et al.|2014| . 


TABLE 2 

Summary of planet and stellar temperatures, goodness-of-fit, and C/O constraints. 


Planet 

HD 209458b 

WASP-19b 

WASP-12b 

HAT-P-lb 

XO-lb 

HD 189733b 

WAS P-17b 

WASP-43b 

T Eq [K] for A b = 1 

1409 

2013 

2517 

1271 

1176 

1169 

1508 

1340 

Aeff M 

6065 

5500 

6300 

5975 

5750 

5040 

6550 

4400 

X^in 

30.8 

1.4 

4.0 

26.8 

35.5 

17.5 

14.4 

24.8 

xlTjN 

1.07 

0.37 

0.58 

0.98 

1.13 

0.81 

0.89 

1.09 

C/O (95% limit) 

< 0.86 

< 0.82 

< 0.80 

< 0.77 

< 0.87 

< 0.92 

< 0.89 

< 0.87 

C/O (99.7% limit) 

< 0.88 

< 0.85 

< 0.87 

< 0.84 

< 0.98 

< 1.20 

— 

— 

( A U) min for C/C»l 

55.9 

13.3 

36.0 

22.9 

12.6 

7.1 

10.0 

6.9 

^max/ ^C/0>1 

1.4- 10 12 : 1 

788 : 1 

7.0 • 10 7 : 1 

99,720 : 1 

553 : 1 

35 : 1 

153 : 1 

32 : 1 


fully accounted for by marginalizing over t he wide range 
of cloud parameters in this model (Section |4.1|4.2|). The 
compositional constraints presented in this section solely 
rely on the HST WFCS observations, and do not depend 
on comparisons of observations take n at different times 
and different instruments. In Section [4~4| I then present 
constraints on the cloud properties using both the “pa¬ 
rameterized particle size and cloud profile model” as well 
as a simplified “gray cloud deck + Rayleigh haze model”. 
HST STIS and Spitzer observations are included to sup¬ 
plement to conclusions on cloud properties. 


4.1. C/O and Metallicity 

The main result of this study is that all eight hot 
Jupiters with detectable near-infrared absorption fea¬ 
tures show a strict upper limit on the C/O ratio at ap¬ 
proximate 0.9. Carbon-rich atmospheric compositions 


(C/0>1) are firmly ruled out for HD 209458b, WASP- 
19b, WASP-12b, HAT-P-lb, and XO-lb — virtually 
no posterior probability exists for atmospheric compo¬ 
sitions with C/O>0.9 (Figure [5b. The remaining two 
hot Jupiters (HD 189733b, WASP-17b, and WASP-43b) 
show similar posterior distributions; however, a low- 
probability tail (<5% probability) remains towards high 
C/O given currently available data (Figure |5l. Mean¬ 
while, the available observations provide virtually no con¬ 
straint on the atmospheric metallicity for any of the plan¬ 
ets. 

The uppers limit of C/O < 0.9 is a robust finding. 
The nested sampling algorithm explores every corner 
of the multi-dimensional parameter space, ranging from 
C/O < 10 -5 (10 5 times more oxygen) to C/O < 10 5 (10 5 
times more carbon), and there is no atmospheric scenario 
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Fig. 7.— Model transmission spectra compared to the published HST WFC3 transit depth measurements (black). Solid lines show the 
transmission spectra for a fiducial clear solar composition atmosphere (yellow), the overall best-fitting model (red), the best-fitting solar 
composition atmosphere wit h clouds (blue), and the best fitting C/0>1 model with clouds (purple ). Plotted obser vations and their 1 — a 
uncertainties a re taken fromp Deming et al. (2013) for HD 20945 8b and XO-lbjHuit son et ahl (120131) for WASP-19 b, Wakeford et al.|(2013) 
for HAT-P-lb, McCullough et al. (2014|) for HD 189733b, |Kreidberg et al.| (|2015|) tor WASF-12b, |Mandell et al.| ([ 20 I 3 ]) tor WASF-lYb, 
and Kreidberg et al. (2014b) tor WAbH-43b. The model spectra tiave a vast number of molecular lines; tor clarity the spectrum has been 
Gaussian smoothed. 
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Fig. 8.— x 2 of all atmospheric models in the nested sampling re¬ 
sults for HD 209458b. The upper limit C/O < 0.9 is extremely 
robust. The nested sampling algorithm explored all corners of 
the multi-dimensional parameter space, ranging from C/O < 10 —5 
(100,000x more oxygen) to C/O < 10 5 (100,000 times more car¬ 
bon). All atmospheric scenario with C/O > 1 are excluded at 
Ay 2 > 55.9, corresponding to a likelihood ratio of 1.4 • 10 12 tol. 
The fit to the data is excellent (% 2 /N = 1.1) C/O < 0.9 and drops 
off sharply near C/O = 0.9 to % 2 /N > 4. 

with C/O>0.9 that results in a good fit to the data. As 
an example, the fit to the data degrades sharply from 
X 2 /N = 1.1 to worse than x 2 /7V = 4 for C to O ratios 


above 0.9. (Figure [ 8 ]). 

The chemical explanation for the strict upper limit on 
the C/O ratio is that the observed water absorption fea¬ 
ture at 1.4 fi m can only exist if substantially more oxygen 
is present in the gas envelope than carbon (Figure [7]). As 
the C/O ratio approaches unity, almost all of the oxy¬ 
gen becomes trapped in CO molecules leaving no oxygen 
left to form H 2 O. The upper limit on the C/O ratio is 
largely independent of the metallicity because the atmo¬ 
spheric chemistry at the high temperatures encountered 
in hot Jupiters sharply changes between C/O = 0.8 and 
C/Oml. 

Methane and water have strongly overlapping absorp¬ 
tion bands in the WFC3 bandpass at at 1.15 and 1.4 /im 
(Figure [9]). Nonetheless, I can unambiguously distinguish 
betweenJJ^O absorption (oxygen-dominated chemistry) 
and CH 4 absorption (carbon-dominated chemistry) be¬ 
cause the relative strengths of the 1.15 and 1.4 /im 
absorption features are vastly different for water and 
methane. Methane absorption bands at 1.15 and 1.4 /im 
are similarly strong and would inevitably result in two 
similarly strong absorption features at 1.15 and 1.4 fi m, 
which is not in agreement with the data. 

Finding a lower limit on the C to O ratio is not pos¬ 
sible because no carbon-bearing species can be inferred 
from the data and a wide range of water abundances 
is consistent with the detected water absorption feature 
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Fig. 9.— Effect of the C/O ratio on the WFC3 transmission 
spectrum of HD 209458b. Colored curves show the best-fitting 
model spectra at C/O=0.59 (blue), C/O=0.75 (red), C/0=1 (yel¬ 
low), C/0=4 (purple), and C/0=17 (green). Best-fitting model 
are obtained by fixing the C/O ratio and optimizing all other model 
parameters including the cloud properties. Water and methane 
have strongly overlapping abs orption bands in the WFC3 band¬ 
pass (Benneke & Seager 2013); the increased transit depth at 
1.15 and 1.4 fi m is due H 2 U absorption for C/O < 1 (oxygen- 
dominated chemistry), but due to CH 4 absorption for C/O > 1 
(carbon-dominated chemistry). All best-fitting models contain a 
thick cloud deck between 2 and 17 mbar to match the relative small 
amplitude of the observed absorption features. C/O < 1 is in good 
agreement with the observations — clouds in the 1-10 mbar regime 
mute the weaker 1.15 jum H 2 O absorption features, while preserv¬ 
ing the stronger 1.4 fi m feature. C/O > 1 is in strong disagreement 
with the data; however, because the methane absorption bands at 
1.15 and 1.4 fi m are similarly strong and would inevitably result in 
two similarly strong absorption features at 1.15 and 1.4 fin 1 , which 
is not observed. 



Fig. 10.— Constraints on the water mole fraction for eight hot 
Jupiters. Colored curves indicate the relative likelihood of the wa¬ 
ter mole fraction at 100 mbar relative to the value expected for so¬ 
lar composition atmospheres (~610 ppm). The water abundances 
of all eight planets are in agreement with solar composition. A 
relatively wide range of water abundances is in agreement with 
the observations due to strong correlation between the water mole 
fraction and the cloud top pressure (see also Figure [T 2 ] ). Account¬ 
ing for clouds in the atmospheric retrieval model, I do not confirm 
the inference of sub-solar water ab undances for HD 209458b and 
HD 189733b (Mad husudhan et al.|2014b| ). 


(Fig. [To]) . The detection of a single absorption feature in 
low-resolution transmission spectroscopy generally pro¬ 
vides weak constraints on the abundance of that absorber 
(Benneke & Seager 2012). The shape and depths of the 
absorption features are predominately determined by the 
scale height and cloud properties, and are only indirectly 
dependent of the absorber abundance. 

4.2. Water Abundance 

The water abundances of all eight hot Jupiters are 
consistent with the value expected for solar com positio n 
(Figures [l 0|. Although previous studies (Madhusud¬ 
han et al. 2014b) have reported low water abundances 
and/or high carbon-to-oxygen ratios on HD 209458b, 
HD 189733b, and WASP-12b, we find that these conclu¬ 
sions are not supported by the data once the full range of 
clouds is considered in the retrieval modeling. It is worth 
noting that the uncertainty in the water abundance is 
predominately driven by the strong correlation between 
water abundance and cloud top pressure (Figures fTol). 
The relatively small depths of the water absorption fea¬ 
tures at 1.4/im as compared to cloud free scenarios clear 
can equally well be explained by scenarios with low water 
abundance and low-altitude clouds (240 mbar) or by so¬ 
lar water abundance scenarios with clouds between 1 and 
100 mbar (Figure[l2|). The correlation between cloud top 
pressure and water abundance explain why lack of clouds 
in the retrieval modeling would lead to the conclusion of 
low water abundance as discussed in the following sec¬ 
tion. 


4.3. Importance of Considering A Wide Range of 
Clouds 


The constraints on the atmospheric composition de¬ 
pend sensitively on the range of clouds and hazes con¬ 
sidered in the retrieval model. It is therefore impor¬ 
tant to account for all plausible cloud scenarios when 
retrieving the atmospheric composition. As an exam¬ 
ple, Figure [Tl] compares the C to O ratio constraints for 
HD 209458b for three different clouds models. With¬ 
out the presence of clouds, the only explanation for the 
relative shallow water absorption feature in the WFC3 
observations of HD 209458b (Fig. [7k is a sub-solar wa¬ 
ter abundance and high C/O ratio (Figure [IT] left). The 
presence of clouds, however, presents an alternative ex¬ 
planation for the observed shallow absorption feature, 
resulting in much weaker constraints on C/O (Figure [Tl] 
center and right). The sensitivity of the atmospheric 
composition constraints is a result of strong degeneracies 
between the atmospheric composition and cloud proper¬ 
ties in low resolution transmission spectra (Figure [12]) . 

4.4. Cloud and Haze Properties 


The observed transmission spectra indicate the pres¬ 
ence of optically thick cloud decks in the mid- 
atmospheres of HD 209458b and WASP-12b (Section 


4.4.1). The clouds decks are found consistently using 
both the “cloud deck + Rayleigh haze”model as well as 
the “parameterized cloud profile and particle size” model. 
HD 189733b may host a similar g ray cl oud deck if star 
spots affect the spectrum (Section |4 4.2); otherwise thin 
haze are required as previously suggested. No clouds or 
hazes are inferred for WASP-19b, XO-lb, WASP-17b, 
and WASP-43b. 
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Fig. 11.— The constraints on the C/O ratio depend sensitively on the range of clouds considered in the retrieval modeling. The three 
panels illustrate the posterior probability density as function of C/O ratio for HD 209458b assuming no clouds (left), accounting for a 
gray cloud deck and Rayleigh hazes (center), and accounting for wide range of clouds using a parameterized cloud profile and particle 
size distribution (right). The plotted posterior probability distributions are marginalized over the model parameters listed in Table [5] 
Completely ignoring the potential presence of clouds leads one to the conclusion of high C/O ratio (C/O = 0.88 d= 0.04) and low water 
abundance (left). The peak for the retrieval model without clouds (left) is narrow because only atmospheres with a narrow range of water 
abundances corresponding to C/O ~ 0.88 fit the relatively small water feature in the absence of clouds. However, only an upper limit 
at C/O=0.89 can be inferred when clouds are considered in the modeling. For the gray cloud deck + Rayleigh hazes model (center), a 
peak at C/O=0.89 remain, but a long probability tail towards C/O —> 0 exists (center). No lower bound on C/O available. The peak at 
C/O ~ 0.88 also disappears in the parameterized cloud profile and particle size model (right). The C/O results presented in this work use 
this parameterized cloud profile and particle size model (right) to fully account for the wide range of plausible clouds on hot Jupiters. 


4.4.1. Cloud Decks on HD 209458b, WASP-12 b, and 
WASPS lb 

It is instructive to investigate the constraints on cloud 
top pressure and haze opacity from the simple “gray cloud 
deck + Rayleigh hazes” model. This two-parameter cloud 
model simultaneously allows for a gray cloud deck at a 
parameterized cloud top pressure as well as “Rayleigh” 
hazes with parameterized opacity composed of small par¬ 
ticles (r p A). Within this model, the HST WFC3 
observations of HD 209458b and WASP-12b strongly in¬ 
dicate the presence of a thick cloud deck. All scenarios 
without clouds are ruled out at greater 99.7% confidence 
(Figures |T2|). Based on the depth of the observed water 
feature size, I infer that the upper cloud deck is in the 
mid-atmosphere between approximately 200 mbar and 
0.01 mbar for HD 209458b and 30 mbar and 0.01 mbar 
for WASP-12b. Stronger constraints on the cloud top 
pressure are not available because the cloud top pres¬ 
sure is strongly correlated with the water mole fraction 
(Figure [12). For HD 209458b, equally good fits to the 
data areobtained for a wide range of atmospheric sce¬ 
narios ranging from low water abundances (1% solar) 
and low-altitude clouds (200 mbar) to high water abun¬ 
dances (10 x solar) and high-altitude clouds (0.01 mbar). 
The clouds in the atmosphere of WASP-12b are between 
0.01 mbar for 5% solar water abundance and 10 mbar 
for 10 times the solar water abundance (Figure [l2|). The 
inferred clouds top pressure for a given water mole frac¬ 
tion is lower for WASP-12 b because the observed water 
absorption feature is weaker as compared to the cloud 
free model transmission spectrum than for HD 209458b 
(Figure [7]). The presence of clouds can inferred from the 
WFC3 observations alone, but is also strongly supported 
by the HST STIS observations (Figure [15]). 

The negative correlations between water mole fraction 
and cloud top pressures in Figure [T2| can be explained as 
follows. High clouds reduce the depths of near-IR water 
absorption features. In contrast, increasing the water 
abundance strengthens the water absorption features for 


a given cloud top pressure. Different combinations of 
water mole fraction and cloud top pressure, therefore, 
result in the same depths of the 1.4 /im water absorption 
feature and similar good fits to the data. 

Modeling the cloud particles using Mie scattering re¬ 
veals that particles with effective diameters above 1 fi m 
are most probable (Figure [l4| ). Large particles natu¬ 
rally explain the gray nature 01 the observed cloud signa¬ 
ture. Both the 1.15 and 1.4 /i m water absorption features 
are diminished, as expected for a sharp, wavelength- 
independent cutoff of the grazing star light below the 
cloud top pressure. It is interesting to note, however, 
that small particles can similarly appear gray in exo¬ 
planet transmission spectra despite their non-gray scat¬ 
tering properties, if the particle density increases sharply 
near the cloud top. A sudden rise in particle number 
density at a given altitude can lead to the cutoff of all 
grazing light beams below that altitude, largely indepen¬ 
dent of the wavelength. The effect of cloud particles on 
exoplanet transmission spectra depends not only on the 
scattering properties of the particles, but also sensitively 
on the vertical distribution of the particles in the atmo¬ 
sphere. 

One caveat for HD 209458 b is that the 1 fi m data 
point at the red end of the HD209458b HST STIS ob¬ 
servations cannot be fit by any model spectrum. Water, 
TiO, and VO are plausible absorbers at 1 /im, but cannot 
account for the increased transit depth without also in¬ 
creasing the transit depth in the surrounding HST STIS 
and WFC3 transit depth measurements. The increased 
transit depth measurement may either be caused by an 
unknown opacity source or systematic effects potentially 
related to the fringing of the STIS CCD may have af¬ 
fected the measurement. 

Similarly, the Spitzer measurement at 3.6 and 4.5 /im 
cannot be fit to better than 1.5 — a when simultaneously 
fitting the HST STIS , WFC3, and Spitzer data points 
(Figure [15)9). While the 1.5 — a deviation of two points 
does not substantially degrade the overall statistical fit 
to the data; it is worth a discussion because they are the 
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Fig. 12 .— Joint constraints on the water abundance and cloud top pressure for HD 209458b, WASP-12b, WASP-19b, and XO-lb. The 
colored shading indicates the relative likelihood of atmospheric models as a function of water mole fraction and cloud top pressure. Black 
contours mark the 68% (lcr), 95% (2cr), and 99.7% (3a) confidence regions, assuming a log-uniform prior on the water mole fraction. 
Scenarios with largely cloud-free atmosphere are near the bottom of each panel; scenarios with high altitude clouds are at the top. Clouds 
are required on HD209458b and WASP-12b at greater than 99.7% confidence. Water abundance and cloud top pressure are highly correlated, 
preventing precise constraints on water abundance and cloud top pressure individually. For WASP-19b and XO-lb, current observations 
do not distinguish between cloud free scenarios and scenarios with low-altitude clouds. Scenarios with and without clouds are in agreement 
with the data. The lower limit on the cloud top pressure depends on the water abundance. 


only measurements red of 3 /im. Sing et al. (2013) and 


Stevenson et al. (2014j) suggested that the entire UV to 
near-lK spectrum is sloped due to small particle hazes. 


Howev er, the new HST WFC3 data by |Kreidberg et ah 


(2015) rule out this scenario unless they are substan- 
tially offset to the HST STIS data. No plausible model 
exists that would simultaneously fit the reported HST 
STIS , WFC3 , and Spitzer data points to within better 
than 1.5 — a. If water absorption increases the tran¬ 
sit depth to 1.46% at 1.4 /im, then the absorption of the 
same water vapor would result in a transit depth between 
1.44% and 1.5% in the Spitzer bandpasses, independent 
of the cloud or haze properties. Possible explanation are 
that the transit depth measurements are low due to ran- 


dom chance or that th ey are affected by star spots (Mc¬ 


Cullough et al.||20 14), WASP-12 b’s stellar companion 
( |Stevenson et ahp 014), stellar variability, or uncorrected 


instrumental effects. 

WASP-17b and WASP-43b show water absorption fea¬ 
tures weaker than expected for a cloud free solar com¬ 


position atmosphere (Figure [7]), which may indicate the 
presence of clouds similar to HD 209458b. However, the 
data is insufficient to rule out low water abundance sce¬ 
narios without clouds for WASP-17b and WASP-43b. 

The observations of WASP-19b and XO-lb do cur¬ 
rently not indicate the presence of clouds or hazes. In¬ 
stead, the observations provide a strict lower limit on 
the cloud top pressure and a strict upper limit on the 
opacity of small-particle “Rayleigh” hazes in the up¬ 
per atmosphere. For WASP-19b, no cloud deck exists 
above the 10 mbar level for solar water abundance and 
0.1 mbar for 50 x solar (Figure [l2| ). Similar, if small 
particle hazes are present, their opacity does not exceed 
<j _r ? 4 oo = 10 -4 m 2 /kg for a solar metallicity atmospheres 
and , cfr ,400 = 1 m 2 /kg for a 50 x solar metallicity atmo¬ 
spheres. Higher values for haze opacity would mute the 
H2O absorption features in the HST WFC3 bandpass 
and result in a sloped visible and NIR transmission spec¬ 
trum, which is not in agreement with the observations. 
Qualitatively the results are similar for XO-lb (Figure 
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Fig. 13. — Observational constraints on effective particle size of 
the clouds on HD 209458b. The most probably scenarios for the 
clouds on HD 209458b are large particles (>1 yn i). Smaller par¬ 
ticles cannot fully be ruled out, however, because a sharp upper 
cloud deck would mask the strongly wavelength-dependent scatter¬ 
ing pr oper ties of small particles in transmission spectra (see also 
Figurefl4|. The assumed particles are Mie scattering particles with 
the complex refractive indices of MgSiC> 3 . 
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Fig. 14.— Observational constraints on the vertical cloud density 
profiles for MgSiO s grains in the atmospheres of HD 209458b and 
HAT-P-lb. Representative cloud density profiles are plotted for 
scenarios that are within a likelihood ratio of 300:1 (Ay 2 < 11.8) 
of the best fit to the observations. For illustrative purposes, cloud 
profiles are shown for near solar composition atmospheres only 
(0.5<M<5; 0.3<C/0<0.6). The cloud density is represented by 
the ratio between the number of condensed MgSiOg molecules and 
the number of background H 2 molecules. The solar abundance 
number ratio Mg/H =4 • 10 -5 is indicated as a vertical dashed 
line. Colors indicate the effective particle radii for each of the cloud 
profiles. Two distinct types of clouds scenarios are in agreement 
with the observations of HD 209458b: 1) a sharp cloud deck near 
1-10 mbar with a cloud free upper atmospheres or 2) thin clouds 
composed of large particles (>1 /in 1 ) reaching to high altitudes. 
Both scenarios reduce the amplitude of NIR molecular feature 
in a largely wavelength-independent way. In contrast, HAT-P-lb 
must host thin, small-particle (< 0.1 ym) reaching to high altitude 
(>1 mbar) to explain the enhanced opacity at visible wavelength. 
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4.4.2. HD 189733b: Small-Particle Hazes or Gray Cloud 

Deck? 

The reported transmission spectrum of HD 189733b 
shows a prominent slope at short-wavelength with in¬ 
creasing transit depth from the near-IR (1 /im) to the 
ultraviolet (300 nm) (Figure \TE\i). This slope has tradi¬ 
tionally been attributed t o scatte rin g of non-absor bing, 
Rayleigh scattering dus t (Lecave lier des E tangs et al. 
20081 Sing et al] 20 111 |P ont et al.||2013|). Recently, 
McCullough et al.| ( 2014| ) suggested that unocculted 
star spots may be responsible for the apparent slope 
in the transit depth, reinterpreting the spectrum of 



Fig. 15.— Ultraviolet to near-infrared transit depth measure¬ 
ments of HD 209458b and WASP-12b compared to model transmis¬ 
sion spectra. Solid lines show the transmission spectra for a fidu¬ 
cial clear solar composition atmosphere (yellow), the overall best¬ 
fitting model (red), the best-fitting solar composition atmosphere 
with clouds (blue). Black circles indicate the tra nsit depth mea- 
surem ent and their 1—0 un c ertainties taken from Knutson et al. 

H Pemin g et al.l d2013h lEva ns et a l. (2015), and Cro ssfield 
2012b) for HD 209458b, ancSing et al. (2013) andKreidberg 
[HUT for WASP-12b. The observed transmission spectra 
indicate tfie presence of optically thick cloud decks on all three 
planets in the mid-atmospheres of HD 209458b and WASP-12b. 
Only the strong 1.4 ym. water band, the cores of alkali lines, and 
the Rayleigh scattering signature blue of 500 nm are detectable. 

HD 189733b’s with a cloud free atmosphere. In this sec¬ 
tion, I revisit the observed spectrum by (1) assuming that 
the reported measurements of the planet-to-star radius 
ratios are not affected by star spots, and (2) accounting 
for the presence of unocculted star spots. In the presence 
of star spots, I find that HD 189733b may have a cloud 
deck in the mid-atmosphere similar to HD 209458b, al¬ 
beit at higher pressure. If confirmed, the cloud properties 
of HD 189733b and HD 209458b could be more similar 
than previously believed. 

First, assuming that the reported measurements are 
not affected by star spots, I find that small-particle 
“Rayleigh” hazes significantly exceeding the opacity of 
molecular Rayleigh scattering must be present in the up¬ 
per atmosphere of HD 189733b, in agreement with pre¬ 
vious studies. Assuming a Rayleigh-like haze opacity 
(cr oc A -4 ), the reference opacity at 400 nm, 400, is 
between ^ 10 -2 m 2 kg _1 for low metallicity atmospheres 
and ^ 10 m 2 kg -1 for 50 x solar composition atmospheres. 
The inferred haze opacity exceeds the expected opacity 
by molecular Rayleigh scattering by a factor of 10 to 
10,000. A strong correlation between haze opacity and 
atmospheric metallicity arises because different combina¬ 
tions of haze opacity and atmospheric metallicity result 
in identical relative transit depths between the visible 
and infrared spectrum. Increased metallicity increases 
transit depth in the near-IR, but the same relative tran¬ 
sit depth between near-IR and visible wavelengths can be 
recovered if the haze opa city is increased simultaneously 
(Benneke & Seager 2012). 
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Fig. 16.— Ultraviolet to near-infrared transit depth measure¬ 
ments of HD 189733b compared to model transmission spectra. 
Black circles in panel (a) indicate the transit depth measure ments 
directly as reported by |Pont et al.| ( |2013] > and [McCullough et al.| 
(2014). Panel (b) shows ( Rp/R *) based on the same measure¬ 
ments, ^ut_coirectedjpr the presence of star spots according to 
|McCullough et al.|j2014) and assuming the best fitting star spot 
fraction and spot temperature. The star spot corrected measure¬ 
ments reveal that the spectroscopic appearance of HD 189733b can 
be similar to the one of HD 209458b and WASP-12b (Figure [l5| . 
Solid lines show the transmission spectra for a fiducial clear solar 
composition atmosphere (yellow), the best-fitting solar composi¬ 
tion atmosphere with clouds and hazes (blue), and a theoretical 
model with only a cloud deck and high altitude hazes (purple). The 
originally reported transit depth measurements (panel a) strongly 
indicate the presence of high altitude, small-particle hazes to match 
the Rayleigh slope between 0.3 and 1. However, an significant off¬ 
set of the WFC3 data is required to match the data. After star 
spot correction (panel b), the measurement do not favor hazy at¬ 
mospheric scenario, but are in better agreement with either a com¬ 
pletely cloud free atmosphere with Na and K removed or scenarios 
with a thick cloud deck at 1-100 mbar similar to HD 209458b. 
Cloud-free scenarios without Na and K simultaneously match the 
HST STIS and WFC3 with all Spitzer significantly offset. Cloudy 
scenarios matches the Spitzer and HST STIS observations with 
the WFC3 offset. Plausible explanations for the suggested offsets 
are variability in unocculted star spot which would give additional 
support to the star spot hypothesis as an explanation for the short 
wavelength slope. 


Alternatively, assuming that star spots affect the mea¬ 
sured transmission, I correct the apparent transit depth 
measurement using the relation 



= £>• 


1-5 



gh //kTphot _ l\ 

e hf/kT spot — l J 5 
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where D is the measured transit depth, R p /R* is 
the geometric planet-to-star radius ratio, 5 is the star 
spot fraction, T spot is the star spot temperature, T p hot 
is the effective temperature of the star’s p hotosphere, 
and / = c/ A is the observing frequency ( McCullough| 


et al. 2014). The corrected transit depth, (R p /R*)' z , 
is illustrated in Figure [H^b) assuming the best fitting 
spot fraction 5 = 4.6% and spot temperature T spot = 
4400 K. After spot correction, the transmission spec¬ 
trum of HD 189733b shows similar trends to HD 209458b 


and WASP-12b. Similar to HD 209458b and WASP-12b, 
a gray cloud deck may be present on HD 189733b that 
reduces the reduced the depth of the water absorption 
feature at 1.4 pm and mute the wide wings of potassium 
and sodium. Molecular Rayleigh scattering in the high 
atmosphere is visible at wavelength shorter than 500 nm. 

5. SUMMARY & CONCLUSIONS 


I have introduced a novel approach to interpret spec¬ 
troscopic observations of planetary atmospheres that 
combines the statistical robustness and exploratory na¬ 
ture of atmospheric retrieval methods 


,, - ers l 

Seager 


Madhusudhan & Seager 2009 Benneke _ w _ f 

Line et al.|2013b ) with the self-consistency and ability to 
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provide physical insight s of complex atmospheric forward 


models ([Burro ws et al. 1997, 2005; Marley et al. 


2002 


Fortney et al. 20Qg nmnMoigTTtTT. 2011b)7TSIEeJ 


than retrieving only molecular abundances as free param¬ 
eters, the new “SCARLET” approach can provide direct 
insights into the elemental composition of the deep atmo¬ 
sphere based on observations of the upper atmosphere. 

In this work, I applied the new SCARLET framework 
to eight transiting hot Jupiters with detected molecular 
absorption features in their near-infrared transmission 
spectra (HD 209458b, WASP-19b, HAT-P-lb, XO-lb, 
HD 189733b, WASP-12b, WASP-17b, and WASP-43b). 
My main findings are as follows: 

1. The C/O ratios in the deep atmospheres of all eight 
hot Jupiters are consistently below 0.9. This result 
strongly indicates that gas envelopes of most hot 
Jupiter are not carbon-dominated (C/0>1). The 
finding of C/O < 0.9 is robust for HD 209458b, 
WASP-19b, WASP-12b, HAT-P-lb and XO-lb 
(Table [2j. A comprehensive exploration of the 
parameter space reveals that any scenario with 
C/O > 1 is excluded at high significance. All 
carbon-rich scenarios (C/0>1) for HD 209458b, 
for example, are excluded at Ay 2 > 55.9 (like¬ 
lihood ratio of 1.2 • 10 12 to 1). The finding of 
C/O < 0.9 relies solely on HST WFCS transmis¬ 
sion spectroscopy, which has been shown to provide 
repeatab le and reproducible measurements (e.g., 
Kreidberg et al. 2014a). No lower limit is provided 


because all C/O ratios smaller than 0.9 are equally 
plausible for all eight planets. 

2. The water abundances for all eight planets are con¬ 
sistent with the value expected for a solar compo¬ 
sition. For HD 209458b, I find that the relatively 
small depth of the 1.4 /im water absorption fea¬ 
ture as compared to a clear solar composition at¬ 
mosphere is due to the presence of a cloud deck, 
not_ du e to a low water abundance as suggested by 


Madhusudhan et al. (2014b). In general, the con¬ 
straints - ^ the water abundance remain weak and 
are driven by the strong correlations between the 
water abundance and cloud top pressure. 

3. Current observations of transiting hot Jupiters pro¬ 
vide no meaningful constraints on the overall abun¬ 
dance of metals in the atmospheres. All metal- 
licity values in the probed range between 0.1 and 
100 times solar metallicity are in good a greeme nt 
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4. Clouds on hot Jupiters are common over a wide 
range of equilibrium temperatures. The transmis¬ 
sion spectra of HD 209458b (T eq = 1400 K) and 
WASP -12 b (T eq = 2500 K) strongly indicate the 
presence of cloud decks in the mid-atmosphere. 
The cloud top pressure is between 0.01 mbar 
and 200 mbar at 99.7% confidence, strongly cor¬ 
related with the atmospheric water abundance. 
HD 189733b may host a cloud deck similar to 
HD 209458b and WASP -12 b, rather than the sug¬ 
gested Rayleigh hazes, if star spots affect the ob¬ 
served spectrum. 

5. An important lesson from this study is that the 
constraints on the atmospheric composition depend 
sensitively on the range of clouds and hazes consid¬ 
ered in the retrieval model. Ignoring the presence 
of clouds or describing the clouds with a more re¬ 
stricted model may produce overly optimistic con¬ 
straints on the atmospheric compositions. In this 
work, I account for a broad variety of cloud proper¬ 
ties and derive reliable constraints on composition 
by marginalizing over free parameters describing 
the particle size and density profile of the clouds. 

6 . The posterior probability distributions of atmo¬ 
spheric parameters obtained from low signal-to- 
noise exoplanet spectra are highly non-Gaussian. 
Parameters, such as the C/O ratio, may have only 
one-sided constraints or long tails at low probabil¬ 
ity. In this regime, presenting the traditional ±1 —a 
( 68 %) uncertainties can be extremely misleading 
because long tails reaching to values far from the 
reported best-fit may substantially change the in¬ 
terpretation of the results. In this work, I discuss 
99.7% confidence limits whenever possible and pro¬ 
vide likelihood ratios to assess our confidence in 
excluding atmospheric scenarios. 


5.1. Implications for Planet Formation 

The upper limit of 0.9 on the C/O ratio has impor¬ 
tant implications for the formation of hot Jupiters. Pro¬ 
toplanetary disk models suggest that oxygen near the 
disk’s mid-plane freezes out efficiently beyond the water 
ice-line, resulting in o xygen-depleted gas between water 
line and CO lines ( Oberg et al.|| 2011 ). Helling et al. 


(2014) showed that the protoplanetary gas between wa- 
ter and CO ice lines should transition to C/O 1 on 
timescales of ^ 3Myrs set by cosmic-ray induced un¬ 
blocking of O 2 and CO. Assuming that most of the giant 
planet envelope directly accretes from this protoplane¬ 
tary gas, the formation of giant planets with carbon- 
enriched envelopes would be a natural consequence. 

The observed upper limit presented here implies that 
either the gas disk has not fully transitioned to C/O 1 
at the time that the planet undergoes runaway accretion 
or the planet’s gas envelope is subsequently strongly pol¬ 
luted by the late accretion of leftover planetesimals. Al¬ 
ternatively, most of the metals in the gaseous envelope 
may be delivered before the planet undergoes runway ac¬ 
cretion. In van Boekel et ah, (2015, in preparation), we 
find that the core of a giant planet can acquire a gaseous 
envelope well before the onset of runaway gas accretion. 
Once this envelope is sufficiently thick for infalling plan¬ 
etesimals to evaporate before reaching the planet’s core, 
the planet’s core stops growing and any accreted plan¬ 


etesimals contribute to the growth of a metal-rich gas 
envelope. Runaway accretion of metal-poor protoplane¬ 
tary gas will eventually dilute this initial envelope with 
vast amount of H/He; however, in this scenario the metal 
composition would be more representative of the solids in 
the protoplanetary disk than of the gas. A metal compo¬ 
sition representative of the solids in the protoplanetary 
disk would naturally explain a low C/O ratio. 

Finally, the observed upper limit on the C/O ratio rules 
out any scenario in which the planet forms inside the 
water snowline and accretes refractory and volatile ma¬ 
terials with ISM-like carbon abundances (van Boekel et 
ah, 2015, in preparation). Such scenarios would lead to 
C/0>1, which is not in agreement with the findings pre¬ 
sented here. Formation within the water ice line remains 
possible, though, if the refractory and volatile materials 
are oxygen-rich similar to refractory and ices materials 
in the Solar System. 

5.2. Future Work 


Observations with wider spectral coverage or substan¬ 
tially more precision and spectral resolution are needed 
to place stringent constraints on the C/O ratios in hot 
Jupiter gas envelopes. Despite the limitations of the cur¬ 
rent data, I was able to place a stringent upper limit 
of 0.9 on the C/O ratios of studied atmospheres due 
to the presence of detectable water absorption in their 
transmission spectra. In the future, the best way of also 
providing lower limits on the C/O ratios is to also de¬ 
tect the absorption of carbon-bearing species such as CO 
and CH 4 . Transmission spectra, in particular, provide 
a straightforward way of determining the relative abun¬ 
dances of H 2 O, CO, and CH 4 by comparing the rela¬ 
tive transit dept hs in the molecular absor ption bands of 
these molecules (Benneke & Seager 2012 ). In this case, 
the SCARLET model can directly deliver the C/O ratio 
from observations covering near-infrared H 2 O and CO or 
CH 4 absorption bands. 

On the planet formation theory side, it will be essen¬ 
tial to understand whether the elemental compositions 
of giant planet envelopes are predominately set by gas 


2011 


Helling et al. 2014) 


phase accretion (Oberg et al. 
or by accretion of solid and ices m the form of planetesi- 
mals (van Boekel et al., 2015 in preparation). For that, it 
will be crucial to understand the relative time scales for 
accretion of an initial gas envelope and planetesimals as 
well as the initialization of the gas runway accretion. It 
will also be critical to understand from modeling whether 
the elemental composition at the bottom of my model¬ 
ing domain at 1000 bars is representative of the entire 
gas envelope or whether the initial formation of a metal 
rich atmosphere and subsequent runaway gas accretion 
can lead to incompletely mixed, stratified gas envelopes. 

Eventually, exoplanet observations may provide a sta¬ 
tistical sample of precise C/O ratios in giant planet en¬ 
velopes ranging from hot Jupiter to wide separation di¬ 
rectly imaged giant planets. The Juno mission — ex¬ 
pected to arrive in 2016 — will obtain the first re¬ 
liable C/O measurement of for Jupiter. Meanwhile, 
disk observation may provide a better understanding of 
the elemental abundances in the protoplanetary gas and 
dust/ices. Together, observations of protoplanet ary disks 
and evolved planets provide the initial conditions as well 
as outcome of planet formation, providing us with the 
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unique opportunity to develop a consistent picture of the 
formation of giant planets. 
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APPENDIX 


EFFICIENT COMPUTATION OF LINE-BY-LINE RADIATIVE-CONVECTIVE EQUILIBRIUM 

Converging to radiative-convective equilibrium generally requires thousands of radiative transfer calculations of the 
entire atmospheric column; henc e previous models have either approxim ated the T-p profiles based on analytical 
solution for (semi-)gray opacities (G uillotl 2010 ; Par mentier fc G uillot|2 Q14|) or b y assu ming co rrelated o pacities across 
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the entire atmospheric column (Burrows et al.|1997 " Marley et al.| 2 (J(J 2 f Fortney et al.|2008 Morley et al. J2013). Here 


I present a novel numerical technique to speed up repeated line-by-line radiative transfer calculations by orders of 
magnitudes — enabling the efficient computation of self-consistent T-p profiles for any atmospheric composition based 
on line-by-line equivalent radiative transfer. The key element of the new numerical technique is a three-dimensional 
radiative-connectivity array, L, describing the radiative link between any two layers in the atmosphere as a function 
of wavelengths. Once L is computed the T-p profile is computed using only a few hundred of spectral points, with 
virtually no loss in the precision compared to the full high-resolution line-by-line radiative calculations with millions 
of spectral points. 

The basic concept of any radiative-convective model is to iterate the T-p profile until the radiative downward flux 
F ^ due to direct star light and infrared reemission is matched by the upward emission flux in the atmosphere. In 
the process, atmospheric layers with temperature gradients exceeding the adiabatic lapse rate — ^ > Y = are 

declared convective and adjusted to the adiabatic lapse rate. For purely absorbing atmospheres, the exact downward 
and upward fluxes are 


oo oo r 

7 t (n) = j F* t xe~ n/,i d\ + j J ttB\ (r) -^-T[ (t; - r) drdX 


0 0 


OO OO T 

in) = J ttB\ (t*) tI (t* -Ti)d\ + J J nB\ (t) (t, - r) drdX 

0 Or, 


(Al) 


(A2) 


where the optical depth t* is the vertical coordinate, r* is the optical depth of the planet’s surface, B\ (r) is the 
Planck function at wavelength A for the atmospheric temperature at optical depth r, and is the stellar irradiation 

at the top of the atmosphere. The term T[ (t* — r) is the slab transmittance 




r)=2 /' 


-Bi-n 




(A3) 


between and r, where r = cos# is the cosine of the inclination towards the upward normal. Conceptually, the 
two terms in Equations Al are the contribution from direct stellar irradiation (1st t erm ) and the downward infrared 
emission from atmospheric layers above r (2nd term). The upwards flux (Equation |A2| ) is composed of the infrared 
emission from the planetary surface or deep interior ( 1 st term) and the infrared emission from the atmospheric layers 
below Tj. Iterating towards convergence, i.e. F^ (u) = F^ (t*) for all r^, is computationally extre mely cos tly b ecause 
each iteration towards convergence requires the evaluation of the double integrals in Equations |A 1 | and |A 2 [ where 
capturing the high-frequency variations in the molecular line opacities requires millions of spectral points. 

The key concept of the new line-by-line radiative-convective model is to break the wavelength integration into parts 
for which the smoothly varying Planck function, B\{r ), is approximate constant, pulling B\ (r) out of the wavelength 
integrals, and storing the piecewise integrated transmission terms in a three-dimensional radiative-connectivity array, 
L, that can be reused for many iterations till t he c onverged temperature pressure profile is obtained. 

As an example, the second term in Equation |A 2 | in discretized atmospheres is 


f N 

UR (a) = j E ~ D ( t b x ) t \ (a - T i) - T \ (n - T i+i) 


dX, 


(A4) 


o •?=* 

where the upward IR flux at level i has contributions from all layers j = i .. .N below the layer i, and B\ ( Tj ) is the 
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Planck function at the temperature T 7 i n layer j. Repeated computations of Equation A4 can accelerated by orders 
of magnitude by rewriting Equation |A4| as 


(Tj , 72 ) Cijki 


j=i k =1 


where 


Cijk 


J [ T \ (n - Tj) ~ T{ (n - T j+ 1 ) 


d\. 


(A5) 


(A6) 




forms the three-dimensional radiative-connectivity array, C and , is the number of spectral elements saved in C. 
Hundreds of spectral elements suffice to perform a line-by-line equivalent evaluation because B (A, T) vari es sl owl y as 
a function of A and no loss information occurs by precomputing CV^. The other term terms in Equations A1 and A2 
can be rewritten in a similar way. 


Once all elements Cijk are computed, the evaluation of Equations A1 and A2 is trivial and convergence to radiative- 
convective equilibrium is achieved within seconds on a modern computer. Since molecular opacities are temperature 
dependent, C^k needs to be repeated 3-4 times, until the model atmosphere is in radiative-convective equilibrium 
and the opacities correspond to the final temperature profile. Overall, the convergence to radiative equilibrium of 
2-3 orders of magnitude faster performing line-by-line computations at each iterative step toward radiative-convective 
equilibrium. 
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